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Abstract 

Single surface plasmon polaritons are excited using photons generated via spontaneous paramet- 
ric down-conversion. The mean excitation rates, intensity correlations and Fock state populations 
are studied. The observed dependence of the second order coherence in our experiment is consistent 
with a linear uncorrelated Markovian environment in the quantum regime. Our results provide 
important information about the effect of loss for assessing the potential of plasmonic waveguides 
for future nanophotonic circuitry in the quantum regime. 
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Surface plasmon polaritons (SPPs) are highly confined electromagnetic excitations cou- 
pled to electron charge density waves propagating along a metal-dielectric interface. A sig- 
nificant effort is currently being devoted to the study of their unique light-matter properties 
and to their use in optoelectronic devices exhibiting sub- wavelength field confinement l|, |2|. 
Most recently there has been a growing excitement among researchers about the prospects 
for building plasmonic devices that operate faithfully at the quantum level. Indeed, the 
hybrid nature of SPPs and the potential for strong coupling to emitter systems via intense, 
highly confined fields {3, 5] offers new opportunities for the quantum control of light 6, 7|. 
The main hindrance to the use of SPPs in practical devices is, however, their lossy char- 
acter. Still, recent work has shown that SPPs can maintain certain quantum properties 
of their exciting photon field, with the demonstration of assisted transmission of entan- 
gled photons 9|, energy-time entanglement 10], quantum superposition 11], quadrature 



squeezing 



121 ]. wave-particle duality [13J and single plasmon detection [141 ] . These results 



suggest that many principles of quantum optics can be transferred to the field of plasmon- 
ics, enabling novel devices such as single-photon switches to be realized 15[] . Despite recent 
progress in using quantum optical techniques to study plasmonic systems, adapting them to 
realistic structures will require a much more detailed understanding of the quantum proper- 
ties of SPPs when loss is present. This is the central focus of our investigation and represents 
an area so far lacking an in-depth experimental study. Understanding how loss affects the 
quantum properties of SPPs may open up a route toward the realistic design and fabrication 
of nanophotonic plasmonic circuits for quantum information processing. 

In this letter, we characterize the effects of loss on the quantum statistics of heralded 
waveguided SPPs. Using single photons produced by type-I spontaneous parametric down 



conversion (SPDC) [4, |l6j, we excite quanta of leaky SPPs in thin metal 



ic stripe waveg- 



uides, one of the fundamental building blocks for plasmonic circuits. [17H20| We measure 
the second-order quantum coherence function, g^(r), Fock state populations and mean 
excitation count rates for a range of different waveguide lengths. We find that the mean 
excitation rate follows the classical intensity rate as the waveguide length increases, but that 
the second-order quantum coherence remains markedly different from that expected in the 
classical regime and keeps a constant value. The measured dependence in our experiment 



is consistent with a linear uncorrelated Markovian environment in the quantum regime 



19]. 



Our results provide important information about the effect of loss for assessing the real- 
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istic potential of building plasmonic waveguides for nanophotonic circuitry that operates 
faithfully in the quantum regime. 

The structures studied consist of 3 fim wide, 150 nm thick gold stripes fabricated by 
electron-beam lithography on glass coated with 23 nm of indium tin oxide. Input and output 
gratings with a periodicity of 680 nm were then etched into the waveguides using focused 
ion-beam milling. To investigate the effects of loss on SPP excitations, the separation 
between gratings was varied between 5 /xm and 30 /xm in steps of 2.5 /im. A scanning 
electron microscope image of a selection of our waveguides is shown in Fig. QJa). These 



asymmetric waveguides support a number of leaky [20j and bound [l7| quasi-TM guided 
plasmon modes. Here, the grating periodicity (see inset of Fig. [U^a)) was chosen to couple 
effectively to the leaky-modes supported by the structure (those with highest field intensity 
at the gold-air interface). Two such modes exist in the waveguides and the lowest order 
mode, shown in Fig. H^b), possesses by far the lowest losses of the complete set of supported 
modes. Indeed, the leaky-modes of asymmetric waveguides are most often used asplasmonic 
circuitry components due to their low losses as compared to the bound modes. 2| 

Fig. [TJ^c) shows a schematic of the experimental setup used to conduct our investigation. 
To create down-converted photon pairs at 808 nm, a 100 mW CW, A=404 nm laser is 



focused onto a BBO crystal cut for type-I SPDC Phase matching conditions |4j, |16[ lead 
to the photons from a given pair being emitted into antipodal points of a forward-directed 
cone with an opening angle of 6° (Fig. [T^c) shows that the antipodal points chosen are 
in-plane). Polarizing beam splitters (PBSs) are placed in the path of the down-converted 
beams to remove any parasitic light with the incorrect polarization. Each beam is then 
injected into a single mode fiber (SMF). A filter is placed on path A to spectrally select out 
the down-converted photons (see Supplementary Information). One of the fibers is directly 
connected to a silicon avalanche photodiode detector (APD A in Fig. QJc)) which monitors 
the arrival of one photon from a given SPDC pair. The second fiber is used to channel the 
other photon from the pair to the sample in the microscope. In our experiment, APD A 
acts as a trigger device, where the detection of a photon in fiber A 'heralds' the presence of 
a photon in the second fiber B (see Supplementary Information). This system allows for the 
generation of heralded photons at a rate of 10 6 s _1 . After collimation, the polarization of the 
generated single photon is adjusted using a half waveplate (HWP) and then focused onto 
the in-coupling grating of the waveguide probed. At the grating, the generated photons are 
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FIG. 1: Experimental configuration, (a) Scanning electron microscope image of a selection of 
waveguide lengths from 5 //m to 20 /zm. All waveguides have been explored in the quantum regime, 
however a further time-domain analysis of the quantum statistics has been performed for the 7.5 /im 
waveguide highlighted by the dashed red box. Inset: detail of one of the in/out-coupling gratings. 
The scale bars denote 5 /mi. (b) Fundamental SPP mode in our stripe waveguide - electric field 
profile along the cross section of the waveguide, calculated using the Finite Element Method (FEM) 
for an infinitely long waveguide, (c) Schematic of the experimental setup including a single-photon 
source stage, waveguide probing stage and final analysis stage. 



converted into SPPs 19[ due to phase-matching conditions (see Supplementary Information). 
These propagate along the waveguide until they reach the out-coupling grating, at which 
point they are converted back into light. The output is selected by an iris and injected into 
a multimode fiber (MMF) for analysis. The dependence of the output light on the input 
polarization (as shown in Fig. EJ^a)) confirms that the collected light originates exclusively 
from out-coupled SPPs. Finally, the multimode fiber directs the output to a Hanbury Brown 
and Twiss interferometer (BS and detectors B\ and B2) used to measure the second-order 



quantum coherence function, g^ 2 \r). 21| All of the photon detection events are time-tagged 



(Hydraharp 400, PicoQuant GmbH). This allows measurements involving detectors B\ and 
B2 to be conditioned on the detection of a photon at detector A, with appropriate delays, 
ensuring that only correlations arising from the injection of single photons into the waveguide 
are measured. 

We first measure the second-order quantum coherence function, g^(r), as a function of 
time delay r between detectors B\ and B2. This coherence function is a measure of the 
correlation of the intensity of a field at a time t — and at a later time t = r for a fixed 
position. By measuring g^(0) for a given field, we can determine whether or not it is in 
the nonclassical regime (g^(0) < 1). In particular, for number states |n), if g^(0) < 1/2 is 
measured in an experiment, we can be confident that the field is within the single excitation 
regime (see Supplementary Information). 

A beamsplitter is used to symmetrically split the field in mode B into modes B\ and B2. 
In this case one can show that the definition of g^ 2 \r) (cf. Equation 1 in Supplementary 



Information) is equivalent to 22 



23] 



( 2 ) M = (^(0)^(^(^(0)) = n BiB2 m 

{E B MEtM^ B2 {r)Et 2 {r)) N Bl N B2 \At)' 

where E + (t) is the electric field operator, T the averaging (integration) time of the mea- 
surement, N BlB2 is the number of coincidence detections at detectors B\ and B2 within a 
coincidence time window At, and N Bl and N B2 are the number of independent detections 
at detectors B\ and B2 respectively. All detections at B2 are delayed by time r. In our 
experiment, we use both an attenuated laser source (A=785 nm) with Eq. [I] used to calcu- 
late the second-order quantum coherence function and the single-photon source, where all 
measurements are conditioned on the detection of a photon in mode A. Therefore in this 



second case we use the conditional form of Eq. CEJ given by [24| 



£\r) = (2) 

^AB^AB 2 

where Na BiB2 is the number of coincidence detections at detectors A, B 1 and B 2 , with 
detections at B\ and B2 occurring within a coincidence time window At centered on the 
detection at A. Nab 1 is the number of coincidence detections at detectors A and B\ within 
the coincidence time window At and similarly for Nab 2 - Na is the number of independent 
detections at detector A. All measurements are taken over an integration time T, which 
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FIG. 2: Intensity dependence and second-order quantum coherence, (a) Normalized coincidence 
rate coupled out from a 7.5/xm waveguide, dependent on the polarization angle of the beam 
incident on the in-coupling grating (red) and theoretically expected cos 2 (0) dependence (black 

(2) 

dashed line), (b) Conditional second-order quantum coherence function, g c (t), for the down- 
converted light in mode B before the waveguide (black), along with the out-coupled light when 
single photons are injected into a waveguide of length 7.5//m (red). The classical limit is illustrated 
by the blue data points, corresponding to the unconditioned second-order quantum coherence, 
</ 2 )(r), for an attenuated laser injected in the waveguide, (c) g^ 2 \r) for the down-converted light 
in mode B (black) and injected into the waveguide (red). Blue: g^'(r) for the attenuated laser 
injected in the waveguide. The inset shows a magnified region of the graph, omitting the red curve. 

does not appear explicitly in Eq. [2j The value of g c (t) at zero time delay provides us with 
a measure of conditioned single-arm statistics in mode B 25j . Additional information about 



the quality of the field intensity corre 



(2) 

ations can then be obtained by measuring g c (r) over 
a range of different time delays r [24j. We have performed such a time-domain analysis in 
our experiment for the out-coupled light from single-photons injected into a fixed waveguide 
length of 7.5/im. The results from the single-photon source before and after the waveguides 
are shown in Fig. E](b) (resp. black and red curves). The time window At used for the 
measurement is 2 ns, and the integration times are adjusted to obtain reasonable error bars. 
The statistics obtained after the waveguide are identical to those of the source itself. One 
can see that the value of g c (0) is < 0.5 in both cases which demonstrates that we are in 
the single plasmon excitation regime. Note that g c (0) is not identically zero in either case. 
The finite value, however, originates solely from accidental coincidences (see Supplementary 
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Information and Ref. [23]). In Fig. |^c) we show the unconditioned g^ 2 '(r) as a function 
of time-delay for the single photon source only and after a 7.5/xm length waveguide (resp. 
black and red). This plot shows the vital role of the detection of photons in mode A for 
the conditional measurements of the statistics of the out-coupled light from the waveguides. 
Without this 'heralding' of the photons, the statistics of the light arriving at the detectors 
are those of a thermal field 24J. It should be noted that the theoretically expected peak of 
g( 2 \0) = 2 for a thermal field is challenging to observe in quantum optics experiments 28], 
as instead of reaching the value of 2, its height above unity is effectively proportional to the 
ratio of the coherence time of the single-photon source to the response time of the detection, 
which in our experiment is ~ 10~ 5 . 

We now turn our attention to the effect of losses in the single excitation regime. As the 
SPPs propagate along the waveguide, the finite conductivity of the metal results in ohmic 



losses, while radiation into the substrate and surface roughness result in radiative losses [29]. 
For a reasonably smooth waveguide surface and a thick gold layer, ohmic losses are the main 
source of damping. Here we operate far enough from the plasmon wavelength for the free 
electron approximation to hold and we thus expect a linear loss model with uncorrelated 
Markovian noise to be valid. In this context we expect that for number states \n) the 
quantum observables that make up transform the numerator of Eq. [T] as n(n — 1) — > 
if n{n — 1) and the denominator as n — > rj n, where i] is the total loss over the length of the 
waveguide 19J. Thus, for this particular loss model the second-order quantum coherence 
should remain unchanged. Note that at the single excitation level, one could anticipate 
situations in which the damping departs significantly from the classical model, for instance 
because of possible correlations between and within the different damping channels, such 
as excited phonons, background ion-cores, electron gas collisions and interband transition 
processes (involving electron-hole pairs). Indeed, closer to the plasmon wavelength, the 
SPP character becomes more electron-like and such effects may become important. This 
regime remains to be investigated. To explore the quantum statistics in the presence of 
loss we first measured the mean excitation rate over a range of waveguide lengths. To do 
this we measured the counts Nb ± at detector B\ for an attenuated laser at a fixed intensity 
and then the conditional counts Nab x at detectors A and B\ for the single-photon source. 
In both cases, the effect of loss from the beamsplitter in the analysis stage was included 
in the overall detection efficiency, enabling us to disregard the data from detector B 2 . In 
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FIG. 3: Count rate statistics of the light out-coupled from the surface plasmon waveguides, (a) 
N Bl for injected attenuated laser as the length of the waveguide is increased. Dashed black line is a 
exponential fit, yielding a propagation length £ = 8.9 ±1.7 /im. (b) Same for Nab 1 at zero delay for 
injected single photons. Here, £ = 9.8±0.6 /im. (c) g c (0) for single photons injected in waveguides 
of varying lengths (green) and unconditioned g^(0) for a laser injected in the waveguides (blue). 
The black dashed line indicates the value found for the single photon source, and the blue dotted 
line the classical limit. 

Fig. |3](a) we show N Bl against waveguide length for the injected attenuated laser. Fig. Mb) 
presents N ABl at zero delay for injected single photons and as the length of the waveguide 
is increased. The Nab 1 trend matches N Bl , providing evidence that the effect of loss on 
the field of the single SPPs is consistent with the classical exponential behavior. The SPP 
propagation length £, denned as the length at which the intensity (mean photon number (n)) 
decreases to 1/e of its original value, extracted from Fig. [3](a) is I = 8.9 ± 1.7 fim, a value 
similar to £ = 9.8 ±0.6 /im obtained from Fig. [3](b). Both values are in good agreement with 
each other, but smaller than the propagation length expected from Finite Element Method 
(FEM) calculations of 16.7/im, due to imperfections introduced by the fabrication of the 
waveguides. 

In Fig. E](c), we show g c (0) for the out-coupled light for injected single photons as the 

(2) 

length of the waveguide is increased. The value of g c (0) for the down-converted photons 
only is plotted as a dashed black line for reference. One can clearly see that indeed the values 
remain unchanged for the lengths investigated. These results therefore provide evidence for 
the validity of a linear uncorrelated Markovian loss model for SPP damping at the single 
quanta level. This complements well and goes beyond previous studies looking into the 
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FIG. 4: Fock state populations of light out-coupled from the waveguides, (a) For conditioned single 
photons injected into the waveguides, (b) A laser attenuated to give on average one photon at 
the detection stage, \a\ 2 ~ 1. In both, the ideal populations are labeled i and the experimentally 
reconstructed populations labeled e. For the attenuated laser, the ideal populations have been set to 
correspond to a weak coherent state with |a| 2 = 1.2. In both plots the errors are calculated from a 
Monte-Carlo approach, propagating the errors from the measured data through the reconstruction 
algorithm. 



preservation of entanglement via localised plasmons [8j, |9| and nonclassicality via long-range 
surface plasmons [3], where elements of plasmon loss were considered. 

We finish our investigation by probing the population structure of the conditioned SPP 
fields propagating along the waveguides. Note that this is the first time, to the best of our 
knowledge, that such a technique has been applied to a plasmonic quantum system. Here the 
number state, or Fock state populations P n for a given state, represents the probability of the 
state to have n excitations. For a single-photon state |1) we have P\ = 1 and P n — 0, Vn ^ 1. 
On the other hand, for an attenuated laser described by the weak coherent state \a>), we have 



the Poissonian distribution P„ 



M 2 k 



-, where \a\ 2 = (n) is the mean excitation number. 
1, we have the first six populations, P — 0.368, 
0.015 and P5 = 0.003, corresponding to the 



For a mean excitation number of | ct | 2 = 
P x = 0.368, P 2 = 0.184, P 3 = 0.061, P 4 
vacuum state |0) and number states |1), |2), |3), |4) and |5), respectively. 

In Fig. IH(a) we show a tomographic reconstruction of the Fock state populations of the 
field out-coupled from a 7.5/im long waveguide, for single photon excitation. In Fig. H^b) 
we show the reconstructed populations for a laser attenuated to give on average one photon 



at the detection stage, \a\ 2 ~ 1. Details of the tomographic reconstruction method 



30] are 



given in the supplementary information. In Fig. H^b), one can clearly see the populations for 
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the out-coupled light from the attenuated laser are consistent with a weak coherent state with 
mean excitation number of ~ 1. On the other hand, the out-coupled light from the single- 
photon source, shown in Fig. Ufa), displays the strong presence of a single population, P x , 
representative of a single excitation. This analysis complements the investigation performed 
in the previous two sections and confirms that we are exciting single SPPs on the waveguides 
when single photons are injected. 

In this work we have used single photons generated by parametric down conversion to 
excite SPPs in metallic stripe waveguides. By measuring the second order quantum coher- 
ence function and the Fock state populations of the light coupled out of the gratings, 
we demonstrated the ability to excite single SPPs. Moreover, the effect of losses incurred 
during propagation of the single SPPs is consistent with the classical exponential behav- 
ior, and does not change the value of g^ 2 \ providing evidence that a linear uncorrelated 
Markovian loss model is valid for SPP damping at the single quanta level. Our results imply 
that building longer and more complex SPP waveguide structures operating in the quantum 
regime is realistic and opens up the possibility for future studies of new types of functioning 
devices based on quantum plasmonics. 
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SUPPLEMENTARY INFORMATION 



1. Second-order quantum coherence function g^ 2 \r) 

For quantized electromagnetic fields propagating in the x-direction with an arbitrary 
lateral beam profile, and represented by the electric field operator E + (x,t), we have at a 
fixed position, x — 0, the following definition [lj 

(E-(Q)E+(Q)) 2 

Here (X) represents the expectation value of the operator X with respect to the initial state 
of the field, i.e. an averaging over ensembles. The average of the intensity of the field is 
assumed to be constant over time, (E~(t)E + (t)) = (E~(0)E + (0)). Throughout we will 
suppress the position dependence of E + (x,t), as x is fixed at zero. At zero time delay, 
t = 0, for n-excitation states \n), we have that (n\ E~(0)E~(0)E + (0)E + (0) \n) = n(n — 1) 
and (n\ E~ (0) E + (0) \n) = n, leading to the relation g^ 2 \0) = 1 — 1/n. In particular, for 
n = 1 (single excitations), we have (/^(O) = 0. Similarly, for n = 2, g^ 2 \0) = 0.5: a 
measured value of ^^(O) between and 0.5 is a confirmation that we are dealing with single 
excitations. On the other hand, for attenuated laser light described by a weak coherent 
state |a) = Xl^Lo e_ ' a ' 2 '^F l n )' wnere \ a \ 2 = ( n ) is the mean excitation number, we have 
(/(^(O) = 1. Moreover, it can be shown using the Cauchy-Schwartz inequality that for 
any classical electromagnetic field, due to the absence of operators and their commutation 
relations for the classical electric field E + (x,t), the numerator in Eq. [1] factorizes to give 
the inequality g^ 2 \0) > 1. Thus by measuring g^ 2 \0) for a given field, we can determine 
whether or not it is in the nonclassical regime (#^(0) < 1). 



2. Spontaneous parametric down-conversion and heralded single photons 

a) Heralded single photons 

The interaction Hamiltonian for type-I spontaneous parametric down-conversion is given 



by 



2 
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H I = h£ a A a B + H.c. (2) 

Here, £ oc X £p> where is the second-order nonlinear susceptibility of the BBO crystal 
in our setup and £ p is the amplitude of the classical coherent laser pump field. In addition, 
a A (a B ) is a creation operator for a photon in mode A (B) and H.c. represents the Hermitian 
conjugate. Taking the initial state of modes A and B to be the vacuum [-0(0)) = |0) AB = 
|0) A |0) B and evolving it according to the Schrodinger equation as 

\m) = e- it6i/h \m), (3) 

we obtain, up to first order in time, the state 

\m) = (1 - k72) \0) a \0) b - i/i \\) A |1) B , (4) 

where fi = £t. By detecting a photon in mode A we remove the first (vacuum) term and 
'herald' the presence of a single-photon state \1) B in mode B, up to first order. By tuning 
the pump laser intensity appropriately, higher order terms can be made negligible in mode 
B, even if the detection in mode A is not photon number resolving. Thus, we can use type-I 
SPDC to produce high-quality single photon states, |1), with larger generation rates than 
currently achieved with emitter-type sources, such as quantum dots [4]. 



b ) Characterization of our single photon source 

In the down conversion process, the phase matching conditions are not perfect in the 
experiment. For this reason,the down-converted light is not monochromatic, but presents a 
spectrum with finite width [5|: this effect is observed with our source as shown in Fig. [U(a), 
where we present its spectral properties. In Fig. QJb) the unconditioned coincidence rates 
Rb 1 b 2 = Nb 1 b 2 /T (where T is the integration time) are shown. No correlation between those 
two arms arises. On the other hand, when we plot the conditioned rates Rab x — Nab%/T 
(Fig. Cijc)) it is apparent that there is a very strong correlation between arm B\ and the 
reference arm A at zero delay. The same occurs for Bi- This figure shows as well that the 
configuration used allows for a single photon generation rate of about 10 6 s -1 . 

In addition, the triple coincidence rate Rab x b 2 = ^ab x b 2 IT is shown in Fig.[]Jd). On first 
thought, one would expect to see a value of zero at zero delay, as a coincidence between A 
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and B\ should indicate the presence of a single photon in the system, and thus forbidding any 
simultaneous detection on B 2 . However, there is a peak in Rab 1 b 2 & t zero delay: this is due 
to the fact that we count coincidences within a finite time window At. In our case, At = 2 ns: 
two "clicks" detected within a 2 ns-wide time window are considered as a coincidence, even 
if they are not exactly simultaneous. For this reason, accidental coincidences are measured, 
at a rate determined solely by the count rates on each detector, the integration time T and 
the time window At. One can show that {(J, if R Bl and Rb 2 are the single count rates at B\ 
and B 2 respectively, the accidental coincidence rate at zero delay for three detectors -R aC c(0) 
is: 



The rate of triples observed at zero delay agree very well with the value expected for solely 
accidental coincidences. Additionally, these accidental coincidences lead to a value of g^(0) 
higher than zero. One can show j(| that the offset on g^(0) due to accidental coincidences 
is: 



By using the values observed in Fig. [T] and the count rates at each detector, one finds 
the value of g^{0) = 0.23 for our single photon source for At = 2 ns, which is the value 
observed in Fig. 2(b) in the main text: the non-zero value of </ 2 ^(0) is solely due to accidental 
coincidences. 



-Racc(O) — AtR ABl R B2 + AtR AB2 R Bl . 



(5) 




(6) 
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FIG. 1: Characterization of heralded single-photon source, (a) Unconditioned spectral intensity 
function for photons in arms A (red) and B (black). Note that a A = 800 nm bandpass filter, 
with bandwidth AA = 30 nm was placed in arm A to spectrally select only a fraction of the 
down-converted photons. This effectively selects out an equivalent bandwidth in arm B when the 
detections are conditioned on A. (b) (resp (c)) Coincidence rate between B\ and B<i (resp A and 
B\). (d) Triple coincidence rate RabxBv 
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3. Fock state population tomography 



Here we provide details of the tomographic method used to reconstruct the popula- 
tions. We used the technique of Zambra et al. [jj to measure the photon statistics based 
on on/off detection. An arbitrary quantum state can be written in the number basis as 
p = J2nmP nm \ n ) ( m l> where the diagonal elements p nn = p n = (n\p\n) give the photon 
number distribution of the state p. Let rj be the efficiency of a given detector, so that 
7] is the probability for a single photon to be revealed and (1 — 77) is the probability for 
it not to be revealed. Thus, the total probability of the detector not giving a 'click' is 
p(rj) = 5^ n (l — f]) n Pn- Consider now a set of such detectors with different efficiencies, r\ v . 
We then have p(rj u ) = ^ n {l — Vv) n Pn, or more compactly written 



Pv 



(7) 



Here, p v can be obtained experimentally and A vn can be set by artificially changing the 
efficiency of the detector, leaving p n as the unknown parameter. Then, by assuming the 
p n 's are negligible for n > n t , where n t is a truncation number and we have at least iV > 
n t detector efficiencies, Eq. [7] is a LINPOS problem j3] and one can use the expectation 
maximization (EM) algorithm, which converges to the maximum likelihood solution. By 
imposing the physical constraint '^ Jn Pn = 1 we have the iterative solution 

Pn + ) = pV ^2 V ~Al rr OO-n' ( 8 ) 

where pn is the value of p n evaluated at the i-ih iteration, f v are the experimental fre- 
quencies of the 'no-click' events for r] = rj u (whose ideal values are p v ) and p v [{pri' > }] are the 
probabilities p v calculated using the reconstructed distribution {p^n} at the i-th iteration, 
i.e. Pu[{Pn}\ = Y^mi^ ~ Vv) n Pn ■ We start the iterations with the unbiased distribution 
pi? = 1/(1 + n t ) and use the experimental results /„ = n 0tU /n u , where for a given detector 
efficiency r] v , n v is the total number of runs (state preparation and measurement) and no M is 
the number of no-click events for these runs. The EM algorithm is then carried out until the 
changes in the population numbers p n between iterations reduce below a given threshold, e. 
For the field out-coupled from the waveguides for the single photon source, we have n u = Na 
and no lU = Na — Nab 1iU Vx- Here, rj x is a loss scaling factor given by r] x = rja/rjo, which al- 
lows us to consider the tomography being performed on the state that enters the detection 
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and analysis stage in our setup (rather than a tomography of the initial state generated, 
as carried out by Zambra et al. {?]]). The loss r/d = 0.55/2 corresponds to detector Bi's 
intrinsic efficiency around the operating wavelength of the field used (A =808nm), combined 
with that of the beamsplitter in front of it. The loss t] Q = Nab 10 /Na is the total loss from 
initial state generation to detection at B\. To measure the coincidences N ab± v we have set 
the coincidence window to At = 2ns. A set of efficiencies are then introduced using an ND 
filter wheel. Here, the efficiencies r\ v = rjdN Bl ,u/NB 10 - For the attenuated laser source, as 
it is not based on conditional measurements at detector A, we set a window of 500 ns for 
a measurement duration every 10 fis and carry out 10,000 runs. Thus n v = 10, 000 and 
n u = 10,000 — Nabi, v i where Nab 1u is the total number of clicks from the 10,000 runs. 
The efficiencies rj u = r] d N BljU /N Bl0 , where r\ d = 0.55/2 is used as before. In Fig. [2] we show 
the dependence of the no-click frequencies f v with the detector efficiencies r\ v measured in 
our experiment for the light out-coupled from the waveguides for the single-photon source 
and that for the attenuated laser. In both plots, background counts were subtracted from 
the singles (Na and N Bl ) and doubles (Nab and Nab^ at the detectors. Using these plots, 
the reconstructed populations from the EM algorithm are found and shown in the main 
text in Fig. 4. To calculate the errors in the populations we used a Monte-Carlo approach, 
propagating the errors from the measured data shown in Fig. [2] through the reconstruction 
algorithm, with Gaussian distributions placed on the values. 
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Iterations Iterations 

FIG. 2: No-click frequencies f v for reconstructing the Fock state populations, (a) No-click fre- 
quencies as detector efficiency rj v is increased for single-photon data, (b) No-click frequencies as 
r/ u is increased for the attenuated laser. In both, the red line corresponds to an ideal single-photon 
and the blue curve corresponds to a weak coherent state with mean photon number \a\ 2 = 1.2. (c) 
and (d) show how the individual populations P n = p n converge as the number of iterations of the 
algorithm increases for single photons and attenuated laser respectively. 
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4. Grating coupling 



Here we provide a theoretical model to describe the transfer of single photons to single 
SPPs via the grating coupling method used in our experiment. At the single-photon level 
only small intensities of the photon field are involved and therefore any nonlinear terms 
in the photon-SPP coupling can be effectively neglected leading to the following linear 
coupling Hamiltonian [9] 

/"OO /"OO 

/ duhutf(u)b(u) (9) 
Jo Jo 

POO 

+ih / du[g(ui)a?(u)b(u) — g*(u)b'(u)a(u)]. 



o 



Here, the a(u;)'s (a^(a;)'s) correspond to annihilation (creation) operators for the photons 
which obey bosonic commutation relations [a (a;), a' (a/)] = 5(u — uj'). Due to the collective 
nature of the electron charge density waves and the frequency regime of our experiment, a 



macroscopic approach for the resulting electromagnetic field is appropriate for the SPPs [10]. 
Upon quantization, they are therefore assumed to behave as bosonic modes. Thus, the 
6(o;)'s (£>' (u;)'s) correspond to annihilation (creation) operators for the SPPs which should, in 
principle, obey bosonic commutation relations [b(uj), tft(u')] = 6(u — u'). In the Hamiltonian 
given in Eq.(9), the first and second terms are the photon and SPP fields' free-energy 
respectively. The last term describes interactions between the two fields, where the coupling 
g{oS) is proportional to the overlap of the field scattered by the grating with the SPP mode on 
the waveguide stripe. Carrying out the overlap integrals for the fields in a particular scenario 
gives the phase-matching conditions for the coupling and the coupling value itself. For 
perfect coupling, i.e. a single photon injected into the grating scatters and couples to a single 
SPP with unit efficiency, we have g(uS) = vr/2 Thus, the geometry of the grating should be 
optimised to achieve a coupling as close as possible to tt/2, with negligible deviation over the 
bandwidth, Aoj, of the incoming photon - in order to avoid significant wavepacket distortion 
and loss during the transfer process. A first approximation for the phase-matching condition 
for photons injected normal to the surface and with TM polarization (with respect to the 
surface plane) is given by |H[ k sp = k g , where k sp is the magnitude of the SPP wavevector 
in the direction of propagation along the waveguide and k g = ^-m is the grating momentum 
(A is the grating period and m is an integer). For the metallic strip in our experiment, 
we use the approximation 



12|, k sp = thi, where e m is the permittivity of gold, 



c 2 l+e m W 2 ' 
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W = 3yU.ni is the width of the waveguide and oj corresponds to a free space wavelength 
of Ao = 808 nm. Taking the m — +1 grating momentum, one finds a grating period of 
A = 802nm. We use this as our starting point and perform FEM simulations to optimise 
the coupling by modifying the height, width, period and number of grooves for the grating 
on the waveguide. We find the optimal period of A = 680 nm. Such a large deviation 
from the approximate result can be explained by the use of deep grooves in our gratings, as 



the phase-matching condition k sp = k g is only approximate for shallow gratings 
perturbations). 
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